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Abstract 

From  the  equations  of  classical  tidal  theory  with 
Newtonian  cooling  (Chapman  &  Lindzen,  1970),  formulae  are 
obtained  for  wind,  temperature  and  pressure  oscillations 
generated  by  thermal,  gravitational  and  lower-boundary 
excitations  of  given  frequency.  The  analysis  is  an 
extension  of  that  of  Butler  &  Small  (1963)  who  formulated 
solutions  of  the  vertical  structure  equation  in  terms  of 
two  independent  solutions  of  the  homogeneous  equation  and 
derived  expressions  for  surface  pressure  oscillations. 

A  comprehensive  formulation  is  presented  for  wind, 


temperature  and  pressure  oscillations  as  a  function  of 
height  with  the  above-mentioned  sources  of  excitation  and 


2 


an  upper-boundary  radiation  condition.  The  formulae 
obtained  are  applied  at  the  surface  leading  to  evaluations 
of  the  surface  oscillation  weighting  function  Wp(z)  which 
weights  the  thermal  excitation  at  height  z  according  to  its 
differential  contribution  to  the  surface  oscillation. 

The  formulae  are  shown  to  simplify  at  heights  above  a 
region  of  excitation  and  evaluations  are  undertaken  of  the 

I 

thermal  response  weighting  function  W^(z)  which  weights 

the  thermal  excitation  at  height  z  according  to  its  ^ 

differential  contribution  to  the  oscillation  at  any  height 

above  the  region  of  thermal  excitation.  Computational 

procedures  are  described  for  obtaining  two  independent 

solutions  of  the  homogeneous  equation  and  results  are 

presented  for  an  adopted  profile  of  atmospheric  scale 

height.  The  problem  of  deriving  the  surface  pressure 

oscillation  due  to  a  tidal  potential  is  briefly  reviewed 

and  results  are  presented  as  an  example  of  the  application 

of  formulae  that  have  been  derived. 
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1.  Introduction 

Since  the  time  of  Laplace  periodic  oscillations  of  the 
atmosphere  as  a  whole  have  provided  a  subject  of  constant 
dynamical  research.  Tidal  components ,  on  account  of  their 
known  periodicities,  have  formed  a  major  part  of  this  study; 
and  reviews  of  investigations  before  1950  have  been  presented 
by  Wilkes  (1949)  and  Chapman  (1951).  Interest  up  to  that 
time  centred  chiefly  around  the  periods  of  free  atmospheric 
oscillation  and  the  so-cailea  resonance  theory,  which 
attempted  unsuccessfully  to  account  for  the  relatively  large 
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magnitude  of  the  solar  semi-diurnal  barometric  oscillation 
in  terms  of  a  free  period  of  very  close  to  12  solar  hours. 
Subsequently,  the  thermal  excitation  of  tides  has  received 
increased  attention  and  was  included  in  the  review  of 
Chapman  &  L indzen  (1970)  from  which  the  equations  of 
classical  tidal  theory  have  been  taken  for  the  developments 
of  this  paper.  In  this  Introduction  some  of  the  main 
features  and  developments  of  classical  tidal  theory  are 
summarized. 

The  theory  of  oscillations  in  a  compressible  atmosphere 
was  initially  developed  as  an  extention  of  that  of  Laplace 
and  others  (Hough,  1S9S)  for  a  liquid  ocean  of  uniform  depth. 
Whereas  for  an  ocean  the  velocity  and  pressure  variations 
were  independent  of  depth,  the  atmospheric  problem  introduced 
a  dependence  on  height  z  which  under  simplifying  assumptions 
could  be  expressed  in  terms  of  the  independent  variable 

X  *  -fnfolpoo  -  |  cl 5/  H(*) 

where  pQ  is  unperturbed  pressure,  pQQ  is  surface  pressure 
and  H  is  atmospheric  scale  height.  The  fundamental 
equation  was 
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ytx 

is  the  ratio  oi'  specific  heats  of  air  and  is 

velocity  divergence. 

The  quantity  h  in  (1.2)  is  a  constant  of  separation 

between  vertically-  and  horizontally-dependent  terms  in  the 

equations  of  motion  and  appears  also  in  the  equation  for  the 

latitudinal  variation  of  the  pressure  perturbation  at  a  given 

height.  This  variation  is  identical  with  that  of  the 

perturbed  depth  of  an  ocean  of  otherwise  uniform  depth  h 

(Laplace's  tidal  equation),  and  it  is  customary  to  refer  to  h 

as  an  equivalent  depth  of  the  atmosphere.  The  choice  of  h 

in  (1.2)  is  limited  to  values  h^’s  which  are  the  eigenvalues 

of  Laplace's  tidal  equation  for  a  given  angular  frequency  of 

oscillation  o’  and  zonal  wave  number  s. 

Use  of  the  alternative  forms  of  wave  representation 

e  and  ,  where  <p  is  longitude  and  t  is 

time,  has  previously  been  discussed  (Groves,  1979) s  the  form 

+  t.) 

that  appears  more  frequently  in  the  literature  is  e  , 

or  being  positive  and  s  =  0,  +1,  +2,  .  .so  that  waves  travel 
westward  if  s  >  0  and  eastward  if  s  <  0,  and  |s|  is  the  number 
of  wavelengths  that  fit  a  circle  of  latitude.  The  form 
ew  ^  will  be  used  in  the  present  paper  (Equ.  2.21). 

The  suffix  n  is  introduced  as  an  identifying  integer  that  is 
assigned  according  to  an  adopted  scheme  of  notation;  then  s', 
s  and  n  specify  a  mode  of  oscillation.  In  the  present  paper 
is  taken  as  fixed  and  only  s  and  n  appear  as  suffices, 


g 

e.g.  hQ.  In  order  tc  reduce  tr.e  number  oi'  suffices  attaching 
tc  various  symbols  the  convention  is  followed  from  the  end  of 
§2  onwards  of  denoting  quantities  dependent  on  s,n  or  summations  of 
such  quantities  by  capital  letters  with  the  suffices  s  and 
n  omitted,  e.g.  q  stands  for  q®:  capital  letters  will  be  used 
exclusively  for  mode-dependent  quantities  and  small  letters 
will  be  used  exclusively  for  mode-independent  quantities. 

One  exception  to  this  rule  will  be  the  continued  use  of  h  for 

s 

hn  and  of  H  for  atmospheric  scale  height. 

The  extension  of  tidal  theory  to  thermal  excitation 
requires  the  introduction  of  a  forcing  term  on  the  right-hand 
side  of  (1.2).  Siebert  (19tol)  presented  the  relevant  equations 
and  investigated  heating  by  water  vapour  absorption  of  solar 

e 

radiation.  If  Jn(x)  denotes  the  modal  rate  of  heating  per 
unit  mass  of  atmosphere,  (1.2)  becomes 


where 


7  -4  ~  a  ■as- 

y_-j  Jn(*> 


0.4.) 


and  g  is  a  constant  acceleration  due  to  gravity. 

For  realistic  profiles  of  scale  height  h,  it  is  necessary 
to  solve  (1.2)  or  (1.5)  numerically.  nefore  1950,  solutions 
of  (1.2)  were  obtained  by  hand  computation  and  differential 
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analysers  for  a  variety  of  scale  height  profiles  in  order  to 
derive  atmospheric  response  curves  (Wilkes,  19^9;  Jacchia  4 
Kopal,  1952).  Integrations  were  carried  out  from  an  upper 
boundary  height  between  125  and  150  km  down  to  the  surface  to 
provide  a  (complex)  solution  y  satisfying  the  upper  boundary 
radiation  condition,  y  was  then  multiplied  by  a  complex  constant 
determined  by  the  amplitude  and  phase  of  the  gravitational 
tidal  potential  through  the  lower  boundary  condition  which 
required  a  vanishing  vertical  velocity  component  at  the  surface. 

Numerical  solutions  of  (1.3)  were  obtained  with  the  aid 
of  an  electronic  computer  by  Butler  &  Small  (1963)  in  an  analysis 
which  showed  heating  by  ozone  absorption  to  be  the  dominant 
generator  of  the  solar  semi-diurnal  barometric  oscillation. 

Butler  &  Small  followed  the  procedure  of  Jacchia  kopal  (1952) 
and  obtained  two  real  solutions  y-^,  yg  of  (1.2)  such 
that  y-^  +  iy^  was  the  required  linear  combination  satisfying 
the  upper  boundary  condition.  The  solution  of  (1.3)  was 
obtained  by  the  method  of  variation  of  parameters  in  the  form 

where 


and  dashes  denote  differentiation.  <*■-* was  determined  as 
previously  to  give  a  vanishing  vertical  component  of  velocity 
at  the  surface  (with  the  tidal  potential  now  disregarded). 


Q  - 


The  availability  of  hign-speed  computers  has  been  a  major 
factor  in  reducing  to  manageable  proportions  the  task  of 
investigating  thermal  atmospneric  tides  for  different  choices 
of  H  and  J®.  Accordingly,  bindzen  (1966)  integrated  (1.3) 
as  it  stands  by  a  method  thac  numerically  applied  the  required 
upper  and  lower  boundary  conditions:  the  corresponding  tidal  wind 
and  temperature  fields  as  wexl  as  surface  pressure  were  also 
computed.. 

The  present  paper  is  concerned  with  developing  the 
analytical  approach  which  has  previously  led  to  (1.5)  and  to 
the  following  equation  for  the  surface  pressure  oscillation 
(Butler  &  Small,  1963) 


X 


o 


o-7) 


On  examining  the  terms  in  (1.7)  it  is  seen  by  (1.4),  (1.6) 

that  the  heating  rate  J®(x)  contributes  differentially  to  the 

surface  pressure  through  the  integral  1^.  Weighting  functions 

Wp(z)  may  therefore  be  defined  for  a  given  mode  which  weight 
s 

JQ(x(z))  in  proportion  to  its  differential  contribution  to  the 
surface  pressure  oscillation  (Groves,  1975).  In  a  similar 
way  thermal  response  weighting  functions  «  (z)  have  been 
introduced  (Groves,  1975*  1976,  1977)  which,  for  a  given  mode, 
weight  J®(x(z))  in  proportional  to  its  differential  contribution 
to  tidal  fields  of  wind,  temperature  or  pressure  at  greater 
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heights.  The  derivation  oi'  a  and  w  has  not  previously 

I  t 

been  given  and  is  included  in  10  and  11. 

The  simplifying  assumptions  under  which  Cl.?)  has  been 
derived  are  detailed  by  Chapman  &  Lindzen  (197C):  the  more 
significant  approximations  are  considered  no  be  the  neglect 
of  (i)  the  Earth's  topography,  (ii)  dissipation  processes, 

(iii)  non-linear  effects  and  (iv)  winds  and  temperature 
gradients  in  the  unperturbed  atmosphere.  The  analysis  developed 
within  this  framework  of  assumptions  has  been  termed  classical 
tidal  theory  (Lindzen,  1968).  The  main  advantage  of  the  theory 
over  more  general  treatments  is  the  considerable  mathematical 
simplification  that  arises  from  the  separation  of  vertical  and 
horizontal  dependences  with  respect  to  eacn  mode  of  oscillation. 
Tidal  motion  in  general  may  then  be  represented  oy  a  summation 
of  such  modes  which  propagate  independently  of  eacr.  other. 

One  notable  extension  of  classical  theory  which  preserves 
separability  has  been  the  inclusion  of  hewtonian  cooling,  i.e. 
a  rate  of  energy  loss  that  is  proportional  to  the  temperature 
perturbation  (Lindzen  &  hchenzie ,  1967;  Lindzen,  196b). 

This  form  of  dissipation  will  be  included  in  tne  present 
analysis. 


Before  proceeding  to  the  treatment  of  vertical  structure* 
series  expansions  of  horizontally-dependent  terms  employed  in 


classical  tidal  theory  will  be  developed.  The  series  coefficients 
introduced  become  the  vertical  functions  involved  in  the  later 
analysis . 

Let  be  an  atmospheric  parameter  which  varies  periodically 

with  constant  angular  frequency  c r  (>0)  and  let 

p.  =  CcrS  9  (.2*') 

where  0  is  colatitude.  We  omit  the  height  dependence  and  write 

Q0(p,4\t  )  =  cohort  +  Cf*,*)  (J-2) 

where  is  longitude,  t  is  time  and  is  a  constant 
having  the  same  physical  dimensions  as  '<Jo.  We  define 

Qc  =  Q*  +  *  Ot 

then 

Qc  “  A}  ^  Qc  e  ^ 

notation:  Superscripts  R  and  I  are  used  to  denote  the  real  and 
imaginary  parts  of  a  complex  number  and  an  asterisk  its  complex 
conjugate. 

Q  is  periodic  in  <i>  in  the  interval  (Q,2i*  )  and  may  be 

V  * 

expanded  as  a  Fourier  series 

OQ 

Q.c  —  21  (  Qc  cc-s  -+-  Qj  -sitYi  C-LS") 

-J=0 

o  p  S3 

where  Qq,  Qg  are  complex.  It  is  convenient  to  replace  Qg 


by  Q-  defined  by 


2  +•  i  Qs 

2G*  -  Qc—  Q5 


0^0) 


c^>  0) 


c^>  0) 


t*.t) 


Then  (2.5)  becomes 


0.  -  Z  Q 


<4 

e 


and  (2.4)  becomes 


o„  -  -CqKi  £  cr*^++'rl'> 


C  2.7) 


^2.8) 


Terms  in  (2.8)  having  s  positive  (negative  )  represent  a  westward 
(eastward)  progression  of  phase.  In  the  special  case  when 
QS  =  0  (s  ^  s  ),  (2.8)  becomes 


Q0  =  JIq  {  Q'4‘R  cnA  t  Q 


where 


t  =  ^of/v  +  t 


O'?) 


(2'0> 


t1  is  then  local  mean  solar  time  if  <r/so  is  equal  to  the  solar 

rate  of  rotation  of  the  Earth  and  t  is  Greenwich  mean  solar  time, 
s 

Q  is  a  function  of  latitude  and  in  classical  tidal  theory 

is  expanded  either  in  terms  of  Hough  functions  ©n  or  of 

that  are  related  to  0s  > 
tyn  n 
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€  gco.V 


-  21 


(2)0 

CQ'0,V)  w.'5) 


where,  for  Q,  =  U,V,W,T  or  P,  the  quantity  is  the  perturbation 
of  the  eastward,  northward,  vertically  upward  components  of 
wind  velocity,  temperature  or  pressure  respectively.  For 
Q=  Qq  refers  to  the  forcing  function  J\  ,  which  is  the 

potential  of  an  applied  force  per  unit  mass  of  atmosphere 
(=  -VJlo),  and  to  JQ,  which  is  the  rate  of  diabatic  heating 
per  unit  mass  of  atmosphere.  The  summation  in  C2.ll)  is 
taken  over  all  members  of  the  set  of  0®  satisfying  Laplace's 
tidal  equation 


J. 

i 

f*  f 

r-F 

Llf-F 

4«X 
9 X 


VP.t3) 


for  -1  <  J-*  ^  1,  where 

j.  =  C3.»4) 

and  ob  is  the  Earth's  sidereal  rate  of  rotation,  a„  its 
o  o 

radius  and  gQ  the  surface  acceleration  due  to  gravity.  0® 
is  the  normalized  eigenfunction  of  (2.13)  cox’responding  to  the 
eigenvalue  h^.  hethods  for  calculating  h®  and  0^  have 
previously  been  reviewed  (Groves,  1979). 

Horizontal  wind  oscillations  depend  on  horizontal  gradients 


of  Sl0  +  PQ/  ^0i  wliere  is  unperturbed  air  density,  and  their 
latitudinal  dependence  takes  the  form  of  (2.12)  where 


o-w*  [  ^  t  il©^ 

Vu* '  f  p  L r  |  J  " 

©  *  =  o-p)H(Wf)p  _  iPlor 

v"  f-p  L  1-0  ^J^"- 


(?-'0 


A  method  for  obtaining  =  U,V)  by  series  expansions 
which  avoid  the  indeterminacy  of  (2.15) »  (2.16)  when  f4.  *  f 
has  previously  been  given  (Groves,  1979). 

Hough  functions  belonging  to  the  same  set  are  orthogonal 


and  when  normalized 


J'0;0>r  =  0 


(n 

(jn.  *  n'  ) 


(2- >7) 


By  (2.4),  (2.7)  and  (2.17),  the  coefficients  in  (2.11)  may 


then  be  obtained  from  Qq  as 


■'2.K 


Qo ( p,4,t )  0JP  e  ol tjfty  (S.I8) 


2**lr 


for  Q  *  W,T,P,J1  and  J.  An  alternative  expression  that 
follows  from  (2.18)  on  replacing  crt  by  <st-^  is 


•iic  _2rr/ar 


c  <rt  ) 


where  by  (2.1) 


=  iQ(Q4  ^o-T- (7.Z0) 

The  height  dependence  of  <^Q  is  expressed  by 

Notation ;  Capital  letters  will  be  used  solely  and  exclusively 
for  quantities  dependent  on  s  and  n,  the  suffices  s  and  n  being 
omitted,  and  for  summations  of  such  quantities.  Capital 
letters  having  suifices  r,j  v.-l,2),  e.g.  A  ,  Crj,  ,  Er, 

Ir,  Kr,  Pr,  Yqj,  Y  .,  and  ZQy  are  also  dependent  on  s 

and  n.  Exceptions  are  made  with  h®,  scale  height  H  and  the  gas 
constant  for  air  R.  . 

By  (2.4),  (2.7),  (2.11)  and  (2.12),  the  expansion  of  is 

Q„  -  4  Rl  0-20 

where  (Q)  denotes  ^  for  Q  =  U,/  or  is  otherwise  omitted. 

The  factors  are  chosen  as  follows 


V 

a4  J 

ip  - 

No/  5f 

i0  = 

L  =  9Jk  “>• 

= 

1  <3  0 

C2.32) 

V 

11 

h 

= 

h  : 

=  j-  Uptitf 

iL  and 


Xjj,  will  be  introduced  in  S§  5  and  4  respectively. 


lb 


3.  Vertical  structure  relations 

The  equations  o  f  classical  tidal  theory  are  taken  from 
Chapman  &  -.indzen  (1970),  where  the  height  dependent  functions 

un’  vn’  wn’  ^Tn’  ^pn*  J^n’  °n  are  relate<i  to  the  quantities 
U,  V,  W,  T,  P,  JL,  J  introduced  in  8  2  by  {2,21)  as  follows 

u-  -  !v  U*  »,  =  5V  V  4W  W 

sT„-  4.  T*  SP..^?*  U-0 

Jl*  J.-'ijJ’ 

A  dependent  variable  Y  is  introduced  here  to  replace  that 
denoted  by  y  in  (1.2),  by  y®  in  il.3)  and  by  yn  in  Chapman 
&  Lindzen  (1970).  We  write 


^  =  («<  -./>’■«  )7  '  0=0 
In  place  of  U,  V,  W,  T,  1  we  begin  by  working  in  terms  of 
N/U  =  -iO 

/  =  V 

-  W-  -ft  C3.3) 

S/T  —  <- U  +  i.Qr)T  +  k  J"  -  l  Xft 

=■  i  [( e*  H/u0)f  ■*-  il] 


where 
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)l  -  JH/  Jz.  (.3-4.^ 

K  -  (*-.)/  V 

in  ( 3 • 3 )  is  a  dimensionless  quantity  related  to  the  rate 
coefficient  of  Newtonian  cooling  denoted  by  a  in  Chapman  & 
Lindzen  (1970).  We  have 

<hr  -  a/<5-V 

Then  Lq,  the  rate  of  decrease  of  the  temperature  perturbation 
Tq,  is  given  by 

L-o/ -£(_  ~  Qr  CM ) 

where  Jl L  is  defined  by  (2.22).  For  a  single  mode  (2.18) 
and  (5.7)  give 

L  *  Qr"^  (3.S) 


Classical  tidal  theory  relates  Y*  to  Y  by 

VQ  -  es\=Qvyj  i«-U,V,vl,T1T)  ^ 

where 

-*v  -rp  z  ^ 

-jLj  ~  U-IA)  - (J-'o) 


JV  -  ( k  4  M(  H'f )  -  V7p 
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Prom  (3.9)  and  (3.10) 


e4x  Y  -  U/H)(Y',/-Yf) 

13..) 

yu *  yV  _  y? 

CS-'O 

yT=  k7tt  c*  +  »Yw 

U>3) 

Y  satisfies  the  vertical  structure  equation 

(Chapman  & 

Lindzen,  1970)  which  may  be  written  as 

-  e4’L[*J/0-aTv 

■  2  c  y  JT J  C3./U) 

where 

3>  --  cf/dx-  -2?  j/dx  -r-fl 

C3.  i5) 

2f  =  i  \ar/(i+iqT) 

(J'O 

R  -  F  - 1  +  r 

C3--7) 

(3.  15) 

-  19  - 


4.  Vertical  energy  flux 

The  rate  of  flow  of  wave  energy  in  a  vertical  direction 
in  a  column  of  constant  cross-section  receives  prominent 
attention  in  tidal  theory  in  connection  with  the  formulation  of 
an  upper  boundary  condition.  As  the  time  average  of  first-order 
flux  terms  is  zero,  products  of  first-order  terms  need  to  be 
retained.  Following  Wilkes  (1949)  the  time-averaged  vertically 
upward  energy  flux  at  a  given  height  in  the  notation  of  the 


present  paper  is 


(A-1') 


where  bracketed  suffices  denote  averaged  quantities.  Hence 


by  (2.4) 


Eh)  ~  a  Jlftly/  CTc  Vic) 


iU.2) 


To  obtain  the  global  average  of  this  quantity  we  first  average 
over  all  longitudes  to  obtain  b,y  (2.7) 

and  then  average  with  respect  to  f*  from  -1  to  +1  to  obtain 
by  (2.11)  and  (2.17) 


(A*-) 


where 


£  *  KICT'W)  o-o 

and  JL^  is  defined  by  (2.22).  By  (3.3)  and  (3*9)  it  follows 
from  (4.3)  that 
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<=  =  Cra0/H;iw[2fr(7>  i5Uli 

=  (<*/«}  Iw.f  7  (JV’/J*  -  i^i4')] 

by  (3.10).  For  x  >  such  that  e~^n  <£  1,  the  term 

in  Si*  may  be  neglected  and  (4.6)  reduces  to 

b  -  o Ji ) i«  t  v  d y7,M  a-7) 

i’he  accuracy  of  this  approximation  is  better  than  1  per  cent  if 
e“^Xfi  a  q#qi,  i.e.  if  height  exceeds  about  70  km.  For  a 
thermal  source  of  excitation,  ^4.7)  is  exact. 

If  in  place  of  (4.1)  we  take 

E(V  ~  CP0  +  f„J*o)(W0  +  JV?o)(t| 

it  can  be  shown  by  (1.1),  (2.4),  (2.7),  (2.11),  (2.17) 
and  (2.22)  that  after  a  short  reduction 

B  =  Rl^fT%  K/He*)jTj(W-;Ji)]  o.<0 

on  putting  H  =*  P0gQ/  pQ-  Hence  by  (3.3),  (3.9)  and  (3.10) 
it  follows  that  (4.7)  holds  without  approximation.  (4.S)  takes 
account  of  the  flux  of  potential  energy  pq  ^0w0  by  replacing 
PQ  by  PQ  +  poJ7o,  where  is  unperturbed  air  density,  and 

the  energy  flux  is  evaluated  with  respect  to  an  equipotential 

* 

surface  for  which  the  vertical  velocity  is  -j70/s0  by  replacing 

9 

WQ  by  W0+J}0/g0.  H  is  then  zero  for  the  equilibrium  tide. 


Let  Yq  -  Yq1,  Y  -j  be  any  two  independent  solutions  of 
the  homogeneous  equation 

-  0 

and  let  Y  =  Y 1  oe  a  particular  integral  of  0.14) ,  then  the 
general  solution  of  (3.14)  is 

1  -  +  V'  isr.z) 

where  ~Q  =  Z- 'zl0>  are  ari)itrar y  constants  and  A  is  defined  by 
notation:  a  a  =  a,  ^  -  o2 

by  definition  the  rtronskian  of  Yq  is 

-UJ0(-0  -  X  A  r*  C.5-  3) 

l*  X. 

and  it  follows  from  (3.14)  and  the  Aoel-Liouv ille  formula  that 

~ .  [X  1 

V  ITO  -  OJ,  (  *<,)  ^'(v)clu  (5'-nX) 

1  L  ■*«  J 

W*  dYo(xo)/dx  may  be  chosen  arbitrarily  as  initial  conditions 

for  the  integration  of  (5.1).  and  nence  by  0.3)  we  may  arrange 

for  w  (x  )  to  be  unity.  Tnis  condition,  will  be  introduced 
o  o  w 

later  (Equ.  12.5)  when  initial  conditions  for  the  numerical 
integration  of  (5*1)  are  considered.  For  the  present  we  note 
that  since  y  is  mode  independent,  wQ(x)  is  also  mode 
independent  by  (5*4). 

By  the  method  of  variation  of  parameters  a  particular 


22 


integral  of  (3. 14)  may  be  obtained  as 

Y'0A,*)=j[  [x M A  X1". 

where  ^ 


S(u)  c/u 


3(u)  =  J"o)  Sn*) 

u5jiu)  =  /['■*•  ^  ar(u)]'u'<-'('u^) 

=  2  i^(u)^  1  V(-u) 


(SiS^ 


(f-0 


and  is  arbitrary.  On  substituting  for  Y1  from  (5*5) 
into  (5-2)  we  obtain 


7bO  =  3>0CO  A  7otx> 


C^7) 


and  hence  that 


where 


dot  °  dx 


Doji*)  =  ~oj<xA)  +  j  X^wSmJu 


^■°l) 


The  suffix  j  (=  1,2)  will  be  used  exclusively  for  the  two 
independent  solutions  of  (3.1) • 

On  substituting  for  Y  in  \3.9)  from  (3*7)  and  noting 
from  (3.10)  that  the  operator  .7^  is  linear  in  d/dx,  we  obtain 
by  (3.8) 
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'iQU)  -  J„(»)A 


O’-  l») 


where  b,y  definition 


From  (5.10)  and  (5«H) 


*=t 


7„j  = 

F 


yU  =  =  7  . 

/Cj  /Oj  Oj 


yt  =  KY?^c*+mv 

'04  «j  °J 

By  (5.9)  we  may  write  (5.7)  as 


O-‘0 


C^'2) 


7(x)  +  j  [7oix)  a  Yotv)j5^d^  -  £>a)a7oi0  ^5l0 


and  (5.10)  as 


-Xg  __  ^ 

VQ(*B)  *  |  [XQ i*sA t 1 i-)]'5t“'a'“ =  -/Vi'A  0*)  ^«0 


on  putting  x  =  xfi.  In  (5.13),  (5.1*0  xA  and  xfi  are  arbitrary 


values . 


I 

J 
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WKBJ  solutions  for  Y  and  E 


WKBJ  solutions  of  (5.1)  are  examined  in  this  section  in 
preparation  for  the  formulation  of  the  upper  boundary  condition 
in  §  7.  By  a  change  of  variable  from  Y  .  to 

7tt.cx)  =  70 -t*.)  C6'0 

J  ^  L  - 


aJ  '  L  Jj  J 

where  £  is  arbitrary,  (5*1)  reduces  by  (3.15)  to 


where 


^  -t-  <Joo  7  to  -  o 

Ax'  J 


Q.  =  F  +  ^  ~ 


Cfc.2^ 


If  G  were  constant,  the  solutions  of  (6.2)  would  be  sinusoidal 
or  exponential.  In  general  G  varies  with  height,  but  if  the 
variation  is  sufficiently  slow,  solutions  of  approximately 
sinusoidal  or  exponential  lorm  may  oe  obtained. 

WKBJ  solutions  may  be  formulated  for  a  range  of  values 
of  x  for  which 

| i  ( £»*/ 0  -  rl  ( -f  I  /  I S I  «  1  (fe  ‘G 

where  suffix  x  denotes  differentiation  with  respect  to  x. 
Solutions  of  (6.2)  are  then  approximated  oy 

[g(x)]  "K  exp  ^  ±  ij  [g(u)]  %  duj  .  Two  independent  wKBJ  solutions 
of  (5.1)  which  satisfy  initial  conditions 


°J  J 

^Vo,(n_  a. 

J  ~  J 

d* 


are  given  by 

[Hjlh-  K  ( £ )]  70jt*)  ”  ^  aj  +  ME)  °i]e»f 

-  [i<*j  +  K.IS^jeapflljIS,*) 

on  using  (6.1),  where 

Ir  (?,x )  ■=  |  Kr  (<0  <W 


x,  -  17  -  i  r. 


r;  -  qk  (-| <«r3r<i)  rs- 


0- 1) 
(t-7  ^ 


The  suffix  r  (  =  1,2)  will  be  used  exclusively  for  terms 
associated  with  the  two  WKBJ  exponential  forms. 

We  define 


fV*'>  - 


l  £t°;L*\  k,u)7.u) 

/ 

-  _ 

K  ,i»)  -  K  (x) 

dLoc  r  °j  J 

r  r 

(6. 10) 


where 

Y  '  ~  3  —  r*  ( t . :  1 ) 

and  Y  . ,  dY  ./dx  are  derived  from  (0.6).  Then  (6.6)  may  be 
0  J  0  Q 

written  by  (6.5)  as 
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V0jlx>>  *  A,j  IS)-e*j»[tIlU,-*')]  +•  ^ij^e->cp[il3l.$^)j  Cfc-12) 

Putting  5  »  x  in  (6.12)  we  obtain 


/\rjlx)+  A  /j  (x)  =  70J  (*)  ^',3> 

From  (6.10)  and  (6.12)  with  x  =  J  '  ,  it  follows  that 


and  from  (6.14)  that 

Ar,(n  __  A,,U) 

Ar,l5')  A,,(.S) 


We  define 


r  i 


L  —°> 

=  __cix _ 

Ara(*> 


0  WX  1 


Crl  -  l 


d  x 


K1  U'X*U) 


then  by  (6.15)  C  ^  is  independent  of  x. 
From  (5.7)  and  (6.12),  we  obtain 


(fc.iv) 


^  -S) 


(6/fc) 


V*/  =  ?,  1,(5,*)]  +  Xl  fe,*)6*/*  J\lz  IS,*)]  C4-'7) 


(x) 

elx 


=  LKt&'fil  If/)  exp  [i  1,1?,/]  4  I K  jx)Y,  e.yp  [fys  *)]  (6.  is) 


where 


IC-II) 


r'r  (5,*)  =  A,  Is) 


Putting  §  =  x  in  (fo.17)  we  obtain 


?ri*,o  -+  ?r,L‘x,-x.)  =  y«.»)  (6.20) 

It  follows  from  (6.14)  and  (6.19)  that 

Trts>>=  Tr(5,>)e=<(>[iIrl'i,S')] 

and  hence  from  (6.17)  that  Y(x)  is  independent  of  $  . 

On  substituting  (6.17),  (6,18)  into  (4.7),  we  obtain  by 
(6.21)  with  5*  =  x 

E  -  E,,  f,  tf'  O  «) 

where 

ErbO=  -Ca./^)K^|Prl5>x)|2ca|i[-5l^f>x)] 

-  -C°o/£)K^)  |PrO,^r 

E'(x)  -  -  *>?.,<.**)}  (W) 

Er  is  an  upward  or  downward  flux  according  to  whether  -K^/h 
is  positive  or  negative.  E1  arises  from  an  interaction 
between  the  two  waveforms  whose  sura  is  Y  in  (6.17). 
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7.  Formulation  of  the  upper  boundary  condition 

Let  xL  be  the  upper  limit  to  the  range  of  x  for  which  a 
solution  for  Y  is  sought.  To  obtain  a  boundary  condition  at 
x  =  x^  consideration  needs  to  be  given  to  the  properties  of 
the  atmosphere  at  x  >  x^.  Following  previous  accounts  (Wilkes, 
19^9;  Butler  &  Small,  1963;  Lindzen,  1966)  we  assume  that 
(1)  the  energy  flux  remains  bounded  as  x-*  °°  and  (2)  the 
radiation  condition  holds,  which  means  that  there  is  no 
incoming  energy  at  large  values  of  x.  The  radiation  condition 
is  usually  applied  on  the  assumption  that  H  and  a,p  are  constant, 
but  the  following  assumption  is  less  restrictive.  We  assume 
that  a  value  of  x^  (  >  xfi)  can  be  found  such  that: 


(i) 

S  =  0  (x  > 

xs) 

IT') 

(ii) 

WKBJ  solutions  are  valid  for  x  ^  x^ 

( iii) 

For  x  5s-  xL, 

either 

(a) 

g1! 

> 

liRl 

17  2) 

or  (b) 

inR 

> 

fl1! 

£.7^) 

By  (5«6)»  (i)  requires  that  fcr  x  Xg 

=  0  C74) 

^A  -  O  L  7-5) 

as  the  two  sources  of  excitation  are  independent.  (7*4)  is 
readily  satisfied  by  taking  Xg  at  a  height  above  the  region 
of  heating  whose  effect  is  being  investigated.  For  example 
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the  combined  effects  of  tropospheric ,  stratospheric  and 
mesospheric  heating  could  be  investigated  o y  taking  Xg  at 
90  km  altitude;  or  the  contribution  of  tropospheric  heating  alone 
could  be  investigated  by  taking  Xg  at  say  15  km  altitude. 

(7*5)  requires,  by  (5.4)  and  (5.6),  that  either  ^  *  0  or 
—  0  to  an  acceptable  order  of  accuracy  for  x  >  Xg. 

Applying  (i)  to  (5*9)  gives 


c*>*s)  (7-0 

Then  (6.19)  with  f  3  xl  and  (7«6)  give 


R^l,?0=  ?  \)  ill) 

Hence  under  (i)  and  (ii),  (6.17)  holds  and  becomes  by  (7.7) 
on  putting  §  »  xL 


+  \  ,  Xs)exP  [l  ^a^c,31)] 

(X.  XL)  (7-s) 


From  (6.8)  we  have 

k,r*  i V 

kr*  -  r?R+  r;r 


k; 


rT-  rK 
>  2. 


<1-  -r1-  rj 


(-7-1) 


I 

Hence  under  (iii)  (a) 

K(x  K*  <  O  (X  C7  '°) 

and  we  choose  rQ  (  =  1  or  2)  s>  ~h  that  <  0  and  K^,  >  0 

o  o 

where  r^  is  given  by  (6.11).  By  (6.7),  (6.25)  and  (7.7)  &r  is 

o 
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unbounded  as  x  -*  unless 

?r  xs)  ^  O  C7-n) 

o 

(7.11)  provides  an  upper  boundary  condition  lor  the  solution 
of  Y  in  the  region  x^  x  4:  x,  . 

O  Xj 

Under  condition  (iii)  (b),  we  have 

K  ,J  K*  >  0  C*  >  x  L )  (7 •  a) 

as  (iii)  (a)  does  not  now  hold;  and  by  (7.9)  that 

( (<*/£)(  !<*/£)  <  O  (x^xj  (7.i5) 

The  requirement  for  a  bounded  energy  flux  excludes  noth 
and  from  being  negative  and  (7.12)  therefore  yields  >>  0 

and  >  0.  By  (7.1$)  the  terms  in  (7.«)  are  then  associ¬ 

ated  with  either  upward  or  downward  finite  energy  fluxes, 
which  by  (6.23)  with  S  =  xl  and  (7.7)  are 

Er=  -  C*o/£)  I 

We  choose  r  (  =  1  or  2)  such  that  /h  )>  0  and  K^,/h  <  0, 

0  ro 

then  by  (7*19-)  is  a  downward  flux  and  E  ,  an  upwai’d  flux. 

o  0 

Under  the  assumption  of  the  radiation  condition,  we  require 
E_  =  0  for  large  x.  Hence  by  (7.14)  F  (x7,x^)  -  0  and  the 

r  xv  Xj  o 

o  o 

upper  boundary  condition  is  again  given  by  (7.11). 

■By  (5.9)  and  (6.19)i  we  can  write  (7.11)  as 
xs  -  _ 

J  *  -y>.)A.Arc»L)  (is) 

*"A  *"  °  c 

On  dividing  by  Ar  2(x^)i  (7.15)  becomes  by  (6.  It-) 


1 


-  '<  i.  - 


5  ^ 

)  fcr  f\  Y  (t.)]  .$£u)dic  ■=  0  0*/O  A 

xA  *-  °  °  J  ° 


C7.'0 


An  upper  boundary  condition  is  provided  by  (7.16)  with  the 
appropriate  choice  of  r  ,  i.e.  rQ=  1  or  2,  such  that  either 

(x)  <  0  (under  (iii)  (a))  or  (x)/h  >  0  (under  (iii)  (b)) 

o  o 

for  x  >xL. 

At  this  stage  we  are  able  to  identify  (1.5)  with  the  general 
solution  (5*15)  subject  to  the  condition  (7.16).  by  (5*9) 


D0j  Lxi_3  =  H0J.  p<a3  +  i  Xj  U 3  5 cm  a-vi  =  ~.OJ  <*  L) 

Hence,  since  S(u)  =  0  for  x  >  ,  (7.16)  becomes 


Cto  -  O 

and,  on  replacing  x  oy  xT ,  (5.15)  becomes 

A  JU 


On) 


a  ■  >«) 


Yc*)  =  ^e(\)A70(x)i  \  A  Y0M]  5^'  A*.  ll  'l) 

°  JC 

Then,  since  C  2  =  ^  anc*-  (7*19)  give 

o 

^  -  SjC,., (,)j  .  [  [lw»^ |u>] -SM  Ju  (7.  JO) 

by  (5.2)  we  may  identify  the  complex  conjugate  of  (7*20)  with 
(1.5)  by  writing 


“t  -  i  =  0 


_  * 

* — * 
<=>0. 


(<x-e  Lf>) 


X,  =  (^K^H 


C  7**2) 

CJ.Z3) 


C«c^o/^)  s'/ 


—  L 


(7*22)  is  in  accord  with  (1.4),  (2.22),  (3.1),  (3.2)  and 
(5»3)»  but  (6.16)  imposes  a  condition  on  y^(x),  y^(x)  for 
x  ^  x^,  which  by  (7.23)  becomes 


Hence 


K‘u  *;},)  C7-*») 

ctoc.  *t> 


«d,-tCyz  «  <z.^p  -il  /]  (x^-*L)  C7. 25"^) 


For  an  atmospheric  'top1  (i.e.  x  >,  x^)  having  constant  scale 
height  H  and  no  dissipation,  i.e.  a^,  *  =  0,  we  have  by 

(3.4-),  (3.18)  and  (6.3)  that  G  is  the  real  constant 
K(H/h)-)^.  If  G  <  0,  (7*9)  and  condition  (iii)  (a) 

require  rQ  =  2  and  (7.25)  gives 

ui+lj z  ^  e*r [-c-gf*]  (7-&) 

If  G  >  0,  (7*9)  and  condition  (iii)  (b)  require  rQ  =  1 
and  (7.25)  gives 

x  +  Ljz  ~  **(*  (-^  (7.27) 

(7.26)  and  (7.27)  are  the  forms  of  solution  introduced  at 
high  level  by  Wilkes  (194-9)  for  this  type  of  atmospheric 
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'top'.  In  general  (7.24)  with  x  =  provides  the  upper 
boundary  condition  in  a  lorn,  that  may  be  applied  to  a  numerical 
integration  scheme  (Lindzen,  1^68). 

It  is  now  possible  to  express  the  upward  energy  flux  at  x 
in  terms  of  the  W&BJ  solution  for  Y.  ^hen  (  7,  11 )  holds,  (7»7)  gives 

*Pr0  .  >•)  0  (~>  *L)  17-2S’) 

and  hence  from  (b.2i)  with  Sj  *  x^,  =  x 

?  C*,x)  =  o  (*-  >  *l)  C7-*1) 

ro 

Therefore  by  (to. 20) 

?r.  (*,*')  -  Y(x)  Cx  >  XL) 

o 

and  the  upward  energy  flux  lin  units  oi  Jly)  may  be  written 
by  (6,24)  as 

F  (->0  )K*  (x)  i  Yum  M  H-3') 

where  K^,/h  <  0. 


I 


-  34  - 


8.  The  relations 

Arbitrary  constants  of  integration  ^02  wbich  were 

introduced  in  S  5  have  since  been  retained  in  the  analysis 
and  in  particular  in  the  upper  boundary  condition  (7.16). 

We  now  turn  to  their  elimination  by  introducing  another 
relation  such  as  a  lower  boundary  condition.  We  take  the 
general  form  of  boundary  condition  expressed  by  (5.14)  and 
defer  its  identification  with  the  lower  boundary  until  the 
last  paragraph  of  this  section. 

From  (5.13),  (5*14)  and  (7.16)  on  eliminating  .c  ,  ^x,), 
we  ob'tain 


.X. 

yoo+J  [y0t*)*X{')]$Lu)<u 

Y«W 

7  (x) 

'ox  ' 

-  0 

[X>  (3C0)aX(u^  5l-)du 

#*.) 

sjQ, 

X2t*8) 

Cc 

c 

rcA 

where  ^  •  U,  V,  W,  T  or  P,  r  = 

1  or  d , 

ana  xA ,  are 

arbitrary.  On  putting  =  xA 

,  (8.1) 

expands  as 

C«-0 


Yt*)  -+■  j  [votx)A7oiu)j5(v')  <tu  = 


Nq^a,x) 

XA 

7s 

(»•=) 


where 


Mq  <•*', *  )  =  X  Cx') A  X  ^ ) 

NqC  *'  x)  =  [<X  A  7  t-x)J  /  [  Cr^  A  7  S<  x  ) 


(X.  5\ 


C*'X) 


An  alternative  form  of  (8.2)  uaj  be  obtained  by  putt  in/' 
'A  S’  XB  ~  "A 


xA  =  xQ ,  xB  *  x.  in  (6.1)  and  expanding  as 


7(^f}X[70ix)A7oiv)]5Mctu  = 

+  j  Mqoa  ,  u)  5(*)cl 


X 

(S-5\) 


(8.8)  also  follows  from  (6.2)  on  using  the  identity 


^Na(>A>x) 


,x)Nq  “ 


V"* 


(u) 


(S.t 


As  xA  is  arbitrary  we  may  put-  xA  -  x  in  (6.2)  and 

obtain 

x 

VcxN^  mq ex,x')'/Q(>)  +  HQ(x,x)j  i\Jq(>  ^  ^.7) 

The  introduction  of  an  upper  boundary  condition  therefore 
enables  Y  to  be  related  to  a  single  Y*  whereas  to 

(3.13)  related  Y  to  two  different  . 

As  the  left-hand  side  of  fb.b)  is  independent  cf  x  a? 

/t 

both  of  which  may  be  assigned  arbitrari  iy ,  the  r:y.r.t-i.a: 
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side  may  te  written  in  two  way.-,  to  give 


[yCt"°  +  jx  [A*a) 

■+  M  ,  (*£ .  w )  $M  du  j 


riy  means  of  the  identity 


(tt.S) 


hMv*>rV’A,“)=  NQi’,^i_/V’,.v> 


l*-P 


(a. 8)  may  be  expressed  as 


n/q (■'*,*'> { ^1q^’ w^lu)‘*u 

*A  # 

+  [70aCx^A'/i1(>'ijJI'SNlQ.(-«,“V5''*W“}  =  'V**,*)''  <V) 

*A  (S.  /  o) 

From  (8.8)  or  (8.10),  Y^(x)  may  be  evaluated  for  a  giver, 
source  function  S(u)  provided  (x  )  is  known  for  some 
particular  Q'  =»  U,  V,  d,  I  or  1,  We  choose  =  W  and 
introduce  as  a  lower  boundary  condition  at  x  =  a  Known 
vertical  velocity  of  the  atmosphere  expressed  non-dimensional ly 
by  W(xA). 
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9.  Formulae  for  the  vertical  dependence  of  oscillations 
As  a  consequence  of  the  linearization  employed  in 
classical  tidal  theory  the  results  of  §  8  show  that  oscillations 
of  an  atmospheric  variable  are  independently  related  to  the 
potential  field,  Si ,  the  diabatic  heating,  J,  and  the  vertical 
motion,  W(xA),  at  the  lower  boundary.  For  a  rigid,  horizontal 
lower  boundary  W(x^)  =  0,  but  more  generally  Earth  surface 
tides  contribute  to  W(x^)  and  the  constraint  imposed  on  air 
motions  by  an  undulating  terrain  introduces  new  modes  that 
have  non-zero  values  for  W(xA),  We  tnerefore  express 
^  =  U,  V,  W,  T  and  P  as  the  sum  of  the  three  above-mentioned 
contributions  by 

a  =  -*■  Qj  ^  Qw  c 9-0 

By  (3.3),  (5«6)  and  (8.8)  with  *  W  we  obtain 

=  i  Lj  )  w (xan  t^'2) 

~*A 

*s 

UT(x^ -  - ; LuuAx)j  M^j ^ 

A  +  ‘■J  S^djW-u^{v'  °^u 


sj  vT  n*) 

J  5 
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Ww(x)  -  Lwc^a>x')  W(xa>)  (v-O 

p  ^  ^ 

Wi W  —  ~ c  L-\A/C:*Ay'x)  -  LJ  ^wsi ^ >4 

xA  J 

-t-  i  (?-7) 

~>c 

Wj  ex)  =  ~LW(*A ,*)j  du 

+  J  J"o^du  CH ■  ^ ^ 


[l+  ^a.rt>)j  l^Cx)  =  -C  LTCXA,>’)  W(xA)  (  Q.  7 ) 

r  ^ 

[it^rw]»=  -LrC^x)[^A)~c|  MWJ^(>AiuQ(u)c/Kj 

+  Vx)Jl(x)-c(  jMTi,c-*(u)^(uidu  Crio) 

-'•x 

[l  +  C  cir(.x)j  Tt>>)  =  t  LtCxa>X^  j  Mw:J.(*Alu)  J(x*')du. 


/  \ 

c[kTcx)-( 


x.s 


MTJ.  (.x  tu)  vJru)  du. 


c^»o 
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[e'HC*)/^ ^  (*>  =  -  c  Lp  x  )  W^A ) 

[ex  HmXJ  g  i-*)  =  -  Lp (*A ,  x ) 


^>2) 


5  *" 

~Lj  ~  -^W 


-  ij  '5MPjic».»*>j7(v')  cM 

X 

fs  H (*)/%] £ CO  =  i-  Lp (*A>* )j  (XA ; U)  Jtv  \  du 

*,4 

r  *5 

-<:i  ^PJ^U)^U)  ^U 


C^-M 


•v  t*  re 


^qq'  U,u')  =  Mq  -jq,  u'> 


Q  ,  * ")  (  ■*• ,  >  ")/  a/q  ( > ,  x  ) 


i.i.'5) 


for  q  =  U,  V ,  wf,  T  or  P  and  1  =  J  or  SI.  An  alternative 
set  of  relations  to  (9.2)  to  (9.14)  may  be  derived  iTom 
(8.10)  in  a  similar  manner. 
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10.  Oscillations  at  the  lower  boundary:  surface  oscillation 
weighting  function 

Expressions  i^r  oscillations  at  the  lower  boundary  may 
be  obtained  from  (9.2)  to  (9.14)  on  putting  x  =  x^. 
Alternatively,  such  expressions  may  be  obtained  from  (6.10) 
with  x  =  xA.  The  latter  procedure  is  followed  as  the 
resulting  equations  simplify  more  readily:  use  is  made  of  the 
relations 


Y^(*)a7ouc*)  =  k‘[Y>)a 

which  follow  from  (3.10), 

I  •+  LpC*.*)  = 

~  UT(^,x)  ~  ^  K 

which  follow  from  (5.12). 
and  (9.16) 


>C(-X>]  =  Y0w(^)h  Yjcx) 

~  -  c* )/& 

(10.,) 

(5*3)  and  (5.11);  and  of 

Hi*)  Nw  (*,*■)/£ 

(Jo-*) 

H(>)  N/w  (>,*)/ £ 

U°-  -3  ) 

From  (6.10)  we  obtain  by 

(3.3) 

"  i'L.yi.x a,xa)  v/Cx/i) 

=  Lu Q*fi,xA)SlLTCfi)  -i  i  [e  (XA)/Rj 

Xs 

^A 

4jtxA)  =  1  [e*A 

34  A 

Ww(xA)  = 

Wjt  C*A>  =  =  O 


Oo.4) 

(lo.s-) 

(io.(.) 

t'o.  7) 
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Tw c -  -  l  lt(»a  v/c »Ay  L i  *  i  oo- 1 ) 

T^C*a>  •=  -k  £nw(xAixa}J2(*a) 


*  ^ 

J* 

-r  [l  +  t  aTl*A)] 

Go.  io) 

TjC*a)  ”  ck|  J*(>A)  —  ^e, A  hx^O J 

*J  N/j-l^,-^)  J"(v.)  <iu|  T  Q  ■+  ‘  a  r(^  {] 

> 

0°'0 

\i  <*A)  -  [«■«  e  *A/  i.  H t»A>J  L?cv , ^ )  Wc>A ) 

e/o.-2) 

Tnc^  *  -Cao/<0  [e5*ANw^*A^(*/0 

-v  t  HOA)j  c/u] 

XA 

CIO.  <3N 

X, 

^T^a)  [ao‘UJoC'*AVt.^J  j  N|^.(>Aiu)  J(V)  civ 

('O'.  w+) 

where 
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Equations  (10.6),  (10.11)  and  (10.14)  show  that  Nj(xa,u) 
may  be  interpreted  as  a  function  of  u  that  weights  J(u)  in 
the  interval  (u,  u+du)  according  to  its  differential 
contribution  for  this  interval  to  oscillations  in  horizontal 
velocity,  temperature  or  pressure  at  the  lower  boundary. 

By  (1.1)  it  follows  that 


^  [ic.ib) 

weights  J(x(z))  in  the  interval  (z,  z+dz)  according  to  its 
differential  contribution  for  this  height  interval  to  the 
lower  boundary  oscillations. 

V/p(z)  has  been  evaluated  for  the  migrating  (s  =1)  modes 
of  solar  diurnal  frequency  designated  by  n  *  +1,  +2,  •  .  +6, 
the  corresponding  values  of  h  being  taken  from  Chapman  & 
Lindzen  (1970).  The  atmosphere  is  assumed  to  be 
non-dissipative  having  a^  =  'f  =*  0  and  hence  by  (3.18)  and 
(6.3) 

<5.  “  ■+■  JH/dx.)  H/ft  -  £  0a,7) 

k  is  taken  to  be  2/7  corresponding  to  V  =  7/3.  The 
adopted  profile  of  scale  height  H  is  shown  in  Fig.  1.  For 
negative  n  it  is  found  that  G  is  real  and  negative,  and 
hence  from  (6.8)  that  ,  K ^  are  imaginary. 


From  (6.16) 
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C  .  are  then  real  if  Y  . (x)  are  taken  to  be  real  when  solving 
J  OJ 

(5*1);  and  hence  Wp(z)  is  real  by  (8.4)  and  (10.1b).  If  the 
heating  maximizes  at  the  same  time  at  all  heights,  i.e.  if 
arg  J(x)  is  constant,  it  follows  from  (10.6),  (10.11)  and 
(10#14)  that  the  phases  of  the  lower  boundary  oscillations 
are  in  quadrature  with  it,  i.e.  they  maximize  earlier  or 
later  by  6  h.  Fig.  2  shows  Wp(z)  on  an  arbitrary  scale 
for  the  first  six  negative  migrating  diurnal  modes  (s  =  1, 
n  =  -1,  .  .,  -6).  Such  modes  are  referred  to  as  trapped 
modes  as  the  generation  of  an  oscillation  by  a  region  of 
heating  decays  in  either  vertical  direction  away  from  the  region. 
The  greatest  contribution  to  an  oscillation  at  the  surface 
therefore  arises  from  atmospheric  heating  Closest  to  the 
surface,  and  there  is  an  exponential-like  reauction  in  the 
contribution  to  the  surface  oscillation  with  the  heigr.t  of 
the  heating. 

For  positive  (n  >  0)  diurnal  modes  a  (z)  is  complex 
and  an  oscillation  at  the  surface  can  be  resolved  into 
components  that  are  respectively  in  phase  and  in  quadrature 
with  the  heating.  The  vertical  structure  of  is 
oscillatory  dividing  the  atmosphere  into  positively  and 
negatively  weighted  regions  (Fig.  $).  From  tne  t-Taph  of 
Wp  for  n  =  1  it  is  seen  that  tropospheric  heating  would 
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Fig.  4  Heal  and.  imaginary  parts  of  w  plotted  on  an 
arbitrary  scale  for  solar  semi-diurnal  modes 

U 

with  s  -  2,  n  ■  2,  .  .  .7.  Hey:  - 

-  Wp.  (s,n)  is  shown  on  each  graph.  The 

imaginary  part  of  (2,2)  is  negligibly  small 
compared  with  the  real  part. 
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I 

be  weighted  with  the  same  sign  at  all  heights,  wnereas 
ozone  heating  which  extends  from  about  20  to  bO  km  is  we in:, ted 
by  two  positive  and  two  negative  regions  thereby  reducing  its 
effectiveness  as  a  generator  of  surface  oscillations.  For 
larger  values  of  n,  vertical  wavelengths  snorten  and  the 
effectiveness  of  tropospheric  heating  also  becomes  reduced 
by  cancellations  between  positive  and  negative  regions. 

Fig.  4  shows  Wp  for  the  first  six  migrating  (s  =  2.) 
semi-diurnal  modes.  For  the  leading  mode  (n  =  2)  W  has  the 
same  sign  at  all  heights  and  hence  tropospheric  seating  and 
stratospheric  heating  combine  in  generating  surface 
oscillations:  although  Wp  decreases  with  height,  ozone  heating 
may  be  shown  to  be  the  main  contributor  to  surface  oscillations 
by  adopting  a  typical  profile  of  tropospheric  and  stratospheric 
heating  (Groves,  1975)*  For  increasing  n,  values  of  h 
decrease  and  become  increasingly  oscillatory.  For  n  =  7, 

Wp  is  almost  identical  with  that  for  n  =  1  in  Fig.  J  as 
both  modes  have  nearly  the  same  value  of  h  (i.e.  O.7O0  and 
0.691  km) . 
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1 1 .  use  illations  above  a  region  ol'  excitation:  thermal 
response  weighting  function  W 

Above  a  region  of  excitation,  expressions  for  the 
oscillations  of  atmospheric  variables  simplify.  with 
x  >  x^,  (o.5)  and  (8.8)  reduce  to 

7 O)  =  oAlx)  2c%)  Ci'-O 

7%  -  LqcvO  2(*a)  (Q.uv.W.TP)  ui.i) 

where  by  (3.3)»  (3.6)  and  (9.13) 

=  W(*„) -  i  Jl(V)  -  j  M 

*A 

Tnen,  by  (3*6)  and  (11.2), 


U t*)  -  :\/o)  =  i  Lui*a,*) 

C"-M 

U'-S) 

[j  -r  t  «TcOj  To>  =  -«■  Lt Oa,x) 

jjexHC‘xV°o]^'*‘)*‘  ~L  Lp  C*a  0) 

l"-7) 

We  define 

Wt(z)  = 

(11.8) 

then  it  follows  from  (11.3)  that  W^(z)  weights  o(x(z))  in  the 
interval  (z,  z+dz)  according  to  its  differential  contribution 
for  this  height  interval  to  $(xA)  and  hence  by  (11. 4)  to 


The  reai  part  of  W^. 
an  arbitrary  scale  for 
s  =  1 ,  n  =  1 ,  •  •  •  6 . 
on  each  graph. 


(nqu.  11.8)  plotted  on 
solar  diurnal  modes  with 
Wj  =  0.  (s,n)  is  shown 
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(11.7)  to  the  oscillations  at  any  p;iven  x  .>  x^.  riots  of 
on  a  relative  scale  are  shown  in  Kips.  S  and  o  for  positive 
and  negative  solar  diurnal  modes  with  the  same  atmospheric 
data  as  in  S  10.  W  is  real  for  a  non-dissipative  atmosphere. 

For  positive  (n  >  0)  diurnal  modes,  is  oscillatory 
dividing  the  atmosphere  into  height  intervals  wh^ch 
make  alternate  positive  and  negative  contributions  to 
the  oscillations  at  a  given  height  above  the  region  of 
excitation.  The  greatest  weights  are  given  to  excitations 
at  the  lowest  heights  on  account  of  the  expor.ential-1  ike 
growth  of  amplitudes  on  propagation  into  air  of  decreasing 
density.  Excitation  by  the  region  of  ozone  neating  is 
much  reduced  by  the  cancellation  of  positive  and  negative 
contributions,  the  reduction  being  greater  for  higher  values 
of  n. 

For  negative  n  (Fig.  o),  nodes  are  characteristically 
trapped  ard  above  a  region  of  excitation  the  greatest  relative 
contribution  to  oscillations  at  a  given  height  arises  from 
the  uppermost  levels  of  the  excitation:  hence  increases 
with  height.  For  n  =  -1  the  trapping  character  is  weak 
and  the  effect  of  a  region  of  heating  extends  over  a 
considerable  range  of  heights. 

Fig.  7  shows  for  the  first  six  migrating  (s  =  2) 
semi-diurnal  modes.  For  n  =  2,  W  changes  sign  at  15  km 
and  therefore  tropospheric  heating  and  stratospheric  heating 


k. 


The  real  part  Wt  of  Wt  plotted  on  an  arbitrary  scale 
for  solar  semi-diurnal  modes  with  s  =  2, 
n  *  2,  .  .  7»  =  0.  (s,n)  is  shown  on 


each  graph 


generate  opposing  contributions  to  oscillations  above  oC  Km, 

For  increasing  n,  become  increasingly  oscillatory  and  the 
extent  to  which  the  contributions  from  tropospheric  and 
stratospheric  heating  to  oscillations  above  60  kjd  combine 
or  oppose  one  another  is  dependent  on  the  structure  of  the 
heating  profile.  For  n  =  ?,  W  is  almost  the  same  as  that 
for  n  *  1  in  Fig.  5. 

The  upward  energy  flux  (in  units  of  |  )  at  x  >  xT 
has  been  expressed  by  (7.31)*  By  (11.1)  this  becomes 

Er,  ("O  =  -(a0/*)K*  (*)  |Mw(^)f  C'l.^) 

O  lC  1 

where  K^,/h  <  0. 

12.  Notes  on  computational  procedure 

Let  (OjX^,)  be  the  range  of  x  for  which  computations  are 

undertaken,  then  in  order  to  evaluate  C  x-,  needs  to  exceed 

r  m 
o 

x^.  Initially  is  not  known  but  values  of  x^  corresponding 
to  150  km  have  been  found  to  be  adequate.  Two  independent 
solutions  of  (5.1)  then  need  to  be  obtained  by 

stepwise  integration.  With  equal  step-lengths  an  undesirable 
loss  of  numerical  accuracy  results  if  solutions  are  exponential 
in  form.  The  difficulty  is  readily  overcome  by  a  change  of 
variable  from  YQ^  to  In  Yq^,  but  a  criterion  is  needed  ry 
which  to  determine  in  advance  of  the  integration  process 
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whether  or  not  such  a  form  of  solution  is  to  be  expected. 

A  criterion  which  has  been  found  suitable  for  treating  the 
solution  as  exponential-like  is  that 

(riI|>|r;Rj 

for  all  values  of  x  in  where  is  defined  by  (6.9). 

Case  (i):  If? I  <  J  |  for  at  least  one  value  of  x  in  (0,x^,) 

Y  .  is  obtained  b,y  stepwise  integration  of  (5.1)  in  the 
^  J 

form  of  two  first-order  equations 


*Y.i  -  V, , 
d  oc  'i 

-  -Ryoj  +  Yj 

02- 3  > 

where  Y-,  •  is  an  auxiliary  variable.  The  integration  is 
-**  J 

started  with  initial  conditions 

y0,(o)  =  =  o 

\0to)  »  o  c<%^o)  ^ 

c/X 

where  6  is  arbitrarily  taken  to  be  the  average  value  of  G  in 
(0,xE).  From  (5.3)  and  (12.4) 

-u>o(o)  -  ,  (j‘2  s) 

Case  (ii):  l^1  |  >  |  ^R|  for  all  x  in  (0,x-^) 


We  write 


1 


yoi  -  yoi  <  k  oj-fc) 

where  K  is  a  constant  which  it  will  be  shown  can  be  chosen 
such  that  (12.5)  holds.  An  auxiliary  variable  X^..  is 
introduced  by 

y.  =  (|2, 7 ) 

then  (5«1)  reduces  to  a  Eiccati’s  equation 

4*u  =  -R  +2f  X  •  -X,* 

d  x 

Initial  conditions  for  Yq1,  YqP  are  taxen  to  be 

V0) -  -2*C-£)4 

d  x 

y0,co  )=2'iC-q.)i  d—°-  (°)  = 

dx 

where  G  is  arbitrarily  taken  to  be  the  average  of  J  in  (U,x-^). 
(12.9)  and  (12.10)  are  appropriate  to  solutions  that  initially 
decrease  and  increase  respectively  with  exponential  form  from 
the  same  initial  value  and  satisfy  (12.5).  by  (12.6), 


(12.7)  the  corresponding  conditions  for  XQ^  are 


dx 

(12. 

XM(o)-0  U2S^10)_ 

dx 

02. 

12s) 

on  choosing 

K  -  i  [*(-§)*] 

(.? 

i3> 

02-9) 

(12  .  iC) 
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The  computational  procedure  would  then  be  to  integrate  (12.8) 
followed  by  integration  of  (12.7)  for  each  of  the  initial 
conditions  (12.11),  (12.12)  and  then  Yq1,  Y ^  would  follow 
from  (12.6)  using  (12.13). 

When  (12.8)  was  integrated  with  initial  conditions  (12.11), 
(12.12)  for  the  case  of  G  real  and  negative,  corresponding 
to  a  non-dissipative  atmosphere  and  solar  diurnal  modes  with 
negative  n,  it  was  found  that  the  solution  with  (12.11)  after 
initially  decreasing  changed  to  one  that  increased  with  x. 

The  reason  for  this  behaviour  may  be  seen  by  writing  (12.8) 
as 

4(X,.-v-)  _  (-q)  -  (x,r  f)1 

d  x 

by  (3.17)  and  (6.3).  For  a  non-dissipative  atmosphere, 

^  =  0,  and  with  initial  condition  (12.12),  variations  in  G 
from  G  result  in  a  solution  X^(x)  ^  [-G(x)]^  and 
dX^/clx  «  o.  With  initial  condition  (12.11)  however  the 
same  variations  in  G  from  G  result  in  a  solution 
X-^(x)  flk  -£-G(x)]^  as  there  is  no  corresponding  change  of 
sign  in  dX^/dx  calculated  from  (12.14).  For  sufficiently 
large  x,  it  was  found  that  X  X^  were  approximately 
proportional  to  X^,  X ^  respectively.  In  principle  the 
solutions  obtained  for  Y  ^ ,  Y  £  are  valid,  independent  solutions; 
but  a  serious  difficulty  arises  because  as  x  increases  they 
become  dependent  solutions  to  the  accuracy  of  computation, 
and  in  subsequent  calculations  all  significant  figures  are 
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lost  through  the  differencing  of  very  large  and  nearly  equal 
numbers . 

The  following  procedure  has  enabled  exponentially 
decreasing  solutions  to  be  obtained  without  loss  of'  accuracy. 
The  solution  with  initial  conditions  (12.12)  is  first 
obtained  as  described.  Then  (12.3)  is  integrated  backwards 
from  x  «  Xj.  to  x  =  0  with  the  initial  condition 

y„(*E ) « -  y,*  bo  (,r'^ 

to  obtain  X-^(x);  and  XQl(x)  is  obtained  from  (12.7)  by 
likewise  integrating  backwards  from  x  =  Xg  to  x  =  0  with 
the  initial  condition 

Xo,  ~ 

By  (12.6),  (12.7)  we  obtain  from  (5.5) 

[X0l(*et  +  Ka(-*e)- (>t)]  Cn.q) 

Hence  by  (5*4-),  (12.15)  and  (12.16) 

w6(o)=  2.  x„  t*e)  t  ^UJ 

In  order  that  (12.5)  should  hold,  we  choose 

k  -  j [2'x„(«n]  -  J  V  «•“ 

ny  (12.6)  the  required  solutions  are 
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^  1 

=  { txr  +  J  t  ‘u]}/[*xia(xto|  *  0^.20) 

o 

Computing  accuracy  may  be  checked  in  the  usual  way  by 
changing  the  step-length  or  by  evaluating  the  Wronskian  by 
(5.3)  and  comparing  the  value  with  that  of  exp  ^2  j  y  duj 
(Equ.  5.4-  and  12.5).  Integration  step-lengths  have  been 
chosen  according  to  the  scheme  in  Table  1.  With  =  K  km, 
an  accuracy  of  1  in  10^  has  been  maintained  in  the  computations 
undertaken. 

Table  1.  Choice  of  step-length  according  to  h.  where  j  is 
step-length  selected  for  0.5  4  h~^  <i  1.25  (km)~*^ 


h  <0:  (-h)'1'*  (kjn)“^  O  -  1,  ,  -  2,  >2 

step-length  >**£  faQ. 

h  >  U:  h-^  (kin)"3*  0-0.5,  0.5-1.28,  1.25-2.5,  2.5-5,  >5 

step-length  i  iii  a  l 


At  the  upper  end  of  the  interval  (0,xE),  the  solutions 
Yqi*  Yq2  are  required  to  approximate  to  the  WKdJ  form  of 
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(6.12)  in  order  to  determine  ^  from  (6.1b).  We  define 

o 

~oi  (XV  12'  ,  -V  V 

Z Kl>1 - J~—. - ' - —  c»4*<*e)  ('3-5') 

*'  <r.M  l  «> 

where  Y  . ,  dY  ./dx  are  obtained  by  numerical  integration  as 

described  above  and  may  not  necessarily  be  WKBJ  solutions. 

In  (12.21),  r  =  rQ,  rQ  being  1  or  2  according  to  the  type 

of  upper  boundary  condition,  i.e.  whether  K1  (x-)  <  0  or 

r  & 

^r  ^  °*  the  va-lues  °Y  Z*r  ^(x)  approximate  to  a 

o  o 

constant  value  for  a  range  of  values  of  x  at  the  upper  end 

of  the  interval  (0,xB),  then  this  constant  value  has  been 

taken  to  be  the  quantity  defined  by  (6.16)  in  which  Yq.., 

dY  ./dx  are  WKBJ  solutions, 
o  J 

Calculations  have  been  undertaken  for  a  non-dissipative 
atmosphere,  i.e.  one  for  which  =  y  =0  and  (10.17)  holds. 
Fig.  6  shows  Yq1,  Yq^  at  heights  from  0  to  150  km  for  various 
values  of  h  for  which  the  Case  (i)  integration  procedure 
applies.  As  h  increases  from  1  to  25  km,  Y  , ,  Y  change 
from  an  oscillatory  to  a  mainly  exponential  form. 

In  choosing  rQ  we  need  to  consider  the  following  relations 
which  derive  from  (6.9)  and  (7.9)  on  noting  that  for  a 
non-dissipative  atmosphere  G,  Gx  are  real 


k:=-  <  =  q 


(£>  o) 


(15.22) 
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Fig.  8  Independent  solutions  ,  Yq  ,  of  the  homcr.eneous 

equation  (5«1).  Values  of  h  are  shown  on  each 

set  of  curves.  Case  U)  integration  procedure 
R  I 

applies.  E  , ,  ,  are  the  real  and  imaginary 

ro1  roi 

parts  of  'L  -^(Equ.  l^.t?l)  which  determines  C  . 

ro'L  1o 

on  approximating  to  a  constant  value  with 
increasing  height.  Calculations  are  for  the 
scale  height  of  Fig.  1  and  a  non-dissipative 
atmosphere  having  V  »  7/8- 


* 
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k,R  =  K*  =  0 

K,1  =  (-£)4-K?x/A<?  t«<0) 

<  =  -(. r'tf+GJ^ 


where  by  (10.17) 


02-^) 


('2.2*-) 


For  all  cases  in  Fig.,  8  &(x^)>  0,  where  x^  corresponds  to 
150  km  altitude,  and  hence  by  (12.22)  with  x  =  x^  the  choice 
rQ  ■  1  complies  with  the  upper  boundary  condition  (iii)  (b) 

of  @  7i  i.e.  (x^J/h  >  0,  K^,(x£)/h  <  0.  Also  it  may  be 

°  o  .  II 

noted  that  condition  (iii)  (a)  of  i  ?  does  not  hold  as  >  u 

by  (12.22). 

Values  of  Z  ,  calculated  from  (12.21)  with  r  =1  are 
ro1  0 

plotted  in  Fig.  8  and  above  100  km  approximate  to  constant 

values  which  determine  G ^  with  an  accuracy  of  the  order 
-5  0 

of  10  in  comparison  with  Cr  ^  =  1 .  The  values  obtained 

for  Cr  ^  depend  on  the  arbitrarily  chosen  solutions 

o  *" 

and  have  no  direct  physical  significance.  ^  enters  the 

o 

analysis  through  (7«16)  which  is  valid  for  any  independent 
solutions  YqH  an<^  the  corresponding  value  of  Cr  ^ . 

Fig.  9  shows  YQ^,  Yq2  for  various  values  of  h  for  which 
the  Case  (ii)  integration  procedure  applies  whereby  one 
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exponentially  increasing  solution  and  one  exponentially 
decreasing  solution  are  obtained.  For  h  =  50  km, 

G(xE)  >  0  and  rQ=  1  as  for  the  cases  in  Fig.  8.  For 
other  values  of  h  in  Fig.  9  G(xE)  <  ”[Gx(xE)/4G(xE)]  2<  0, 
and  hence  by  (12.23)  <  0  and  the  appropriate  choice 

of  rQ  is  rQ=  2  under  condition  (iii)  (a)  of  §  7,  i.e. 

K*  <  0  and  K*  >  0. 

For  all  cases  in  Fig.  9,  1  =  0  and  Z^  1  decreases 

o  o 

at  sufficiently  great  heights  to  very  small  values  indicating 

that  cr  p  ~  Cr  1  ^cwever  multiplied  by  Y  „  in  the 

o  o  0<~ 

product  C  i  a  Y  and  as  Y  increases  indefinitely  it  does 

not  necessarily  follow  that  the  product  is  negligibly  small. 

D 

The  orders  of  magnitude  of  Zf  L  at  130  km  are  shown  in  Fig,  9 

oA 

and  when  multiplied  by  Yq^  the  product  is  seen  to  be 
negligibly  small  compared  with  unity.  The  same  result  may 
be  established  analytically  using  the  relation 

vy  y  _  ,  (j 2.  as) 

dx.  d  Ol 

which  follows  from  (5.3),  (5*^)  and  (12.5).  Then  by 
(6.16)  and  (12.25) 

A  X  -  K. ,  7^  -  cl  Yoa/  dxj  (ii..t) 

If  rQ  =  1,  as  for  the  case  h  =  50  km,  G(x)  >  0  for  values 
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oi'  x  close  to  x,,  and  by  (1 2.2a  j  and  (12.26) 

Hj 


I 


Hence  if  Y^,  dY0Vdx  are  lar.^e  compared  with  unity  the 
right-hand  side  of  (12.27)  is  correspondingly  small  compared 
with  unity  as  G  >  0.  If  rQ  =  2  as  for  the  other  values  of 
h  in  Fig.  9,  G(x)  <  0  for  values  of  x  close  to  x^.  and  by 
(12.23)  and  (12.2b) 


|<VX 


Hence  for  large  Yo2,  dYo2/dx  (which  are  positive  for  tne 
exponentially  increasing  solution)  the  right-hand  side  of 
(12.28)  is  correspondingly  small  compared  with  unity  as 
K*  >  0. 


13.  Surface  pressure  oscillation  due  to  a  tidal  potential 

The  variation  of  the  surface  pressure  oscillation  F^(u) 
with  tne  equivalent  depth  h  has  teen  extensively  investigated 
for  different  profiles  of  scale  height  H  in  connection  with 
the  theory  of  atmospheric  resonance.  we  follow  earlier 
procedure  (Wilkes,  194-9;  Gacchia  ft  nopal,  1932;  Giebert, 

19bl ;  nutler  &  Gmail,  19b3;  Giwa,  19bb)  and  tone  the  ratio 


u  _ 


v^)  =  .o  os.>) 

A 

where  P^  is  Riven  by  (10.15)  and  Pg  is  the  dimensionless 
height-dependent  function  for  the  equilibrium  tide. 

For  the  solar  semi-diurnal  oscillation  it  is  now 
recognized  taut  the  effect  of  tnermal  excitations  dominates 
the  gravitational  tide,  but  for  the  lunar  semi-diurnal 
oscillation  this  is  clearly  not  the  case  and  the  evaluation 
of  77  for  the  appropriate  value  of  h  (  =  7.0?  km)  has  seen 
of  continuing  interest.  rsy  direct  numerical  integrations 
of  the  classical  tidaL  equation,  results  for  io  with 
h  =  7.07  km  have  been  shown  to  oe  extremely  sensitive  to  the 
choice  of  basic  atmospheric  properties,  namely  the  profiles 
of  scale  height  and  Newtonian  cooling  i Chapman  h  Lindzen, 

1970).  In  descriptive  terms,  the  sensitivities  arise  from 
the  setting  up  of  a  dependence  on  multiple  reflexions  between 
horizontal  surfaces  at  different  heights  which  are  critically 
dependent  on  basic  atmospheric  properties.  ^hen  the  zonal 
winds  and  latitudinal  temperature  gradients  of  a  realistic 
atmosphere  are  introduced  the  reflecting  surfaces  are  no  longer 
horizontal  and  by  integration  of  the  non-classical  equations 
multiple  reflexions  and  resulting  sensitivities  are  found 
to  oe  largely  removed  (Lindzen  &  Hong,  1974).  we  shall 
nevertheless  develop  the  analysis  for  t7  by  classical  theory 
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in  view  of  tne  earlier  attention  that  it  received,  rut  for 
some  range  of  values  of  h  incluairg  h  =  7.07  *m  it  is 
recognized  that  the  results  are  not  realistic. 

By  definition  the  static  equilibrium  tide  has  zero 
horizontal  components  of  velocity,  and  hence  from  if.  3./  and 
(5.12)  =  0.  From  (5.5)  we  then  obtain 

tf  -  -  e' XJl  a./  H  (15.2) 

By  (10.13)  with  x  •  u  and  i13.2)  with  x  =  0,  we  may  write 
(13-1)  by  (10.13)  ana  <12.y)  as 

ia(rK)  ^  [N/w  ( 0,  0)  f  i|  [0,  Ju  Hl0)/ft  ■*) 

C 

where  the  very  small  dependence  of  Jl(x)  on  x  has  been 
neglected.  by  (6.4),  we  may  men  write  (13. x)  as 

t;  ( r0 ,  a )  =  £h-± tL(_°)  (i3.. 

^roAY0W(0)  ^ 

where 

Vooj  =  V  (o)  +  if  5V.  C'V«r) 

An  infinite  response  therefore  arises  for  values  of  h 
which  satisfy 


cr  *  r<>) 

G 


0 


(13.  U  ) 


Oscillations  lor  which  (13. o)  holds  are  relerred  to  as  tree 
osc il lat ions . 

We  proceed  on  tne  assumption  tnat  the  atmosphere  is 
non-diss  ipative ,  i.e.  that  a  =  '•f  =  0  and  tnat  ilG.l"), 

(1  4.22)  to  (1  2.24;  no  la.  :  .en  oy  (p.o;  =  0,  :\v  (l^.S) 

Yooj  *  Yoj^  and  aay  De  written  Dy  (.3.10/  and  (5.11) 

as 


! 

^  V(0)]h< 0)H 

CV  a 

C 

(IV/'l 


when  Case  (i)  integration  procedure  applies,  we  obtain 
on  introducing  (12.4)  into  (13.7) 


V  ( o ,  ^  ) 


-i 


-  I 


-  U  -C,  5'*)4/Hio) 


(I5.S) 


when  Case  (ii)  integration  procedure  applies  such  that 

Yq^(x)  nas  an  exponentially  decreasing  form  and  r.as 

an  exponentially  increasing  lorm,  we  have  that  ^  =  0. 

1  o1 

Hence  (13.7)  gives  by  (12.6)  and  (12.7) 


0^,-0]"=  I  -  [if  -X„(o/jl/HlO)  0^1) 

In  the  special  case  when  x^  =  0,  i.e.  when  wKHJ 
solutions  hold  for  all  x,  G(x)  which  is  real  by  (10.17) 
has  the  same  sign  at  all  heights  otherwise  (b.4)  would  be 
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invalidated.  Under  Case  (i)  integration  procedure  we  have 
that  G(x)  >  0,  and  from  (12.88)  it  follows  that  the  choice 
r  =  1  complies  with  the  upper  boundary  condition  (iii)  (b) 
of  g  7,  i.e.  kJS  (x)/h  >  0,  K*,(x)/h<  0  (as  h  >  0  by  (10.17; 

in  which  K  +  dH/dz  >  0  for  a  realistic  atmosphere).  (6.16) 
now  holds  for  x  =  0  and  by  (12. 4)  gives 

C„  =  -  C  K2(  o)/  £  z  U'S-to) 

Hence  (13.8)  becomes 

£7(!,n>J'' -  I  -  [3  -i  #Ca(o>j-ft/  H(OJ 

=  <  ~^  +  f|fo))  +  i[Q'(Q)]i)^/H(0)  05  ,l) 

by  (12.22).  Under  Case  (ii)  integration  procedure  we  have 

G  <  0  and  by  (12.23)  it  follows  that  the  choice  rQ  =  2 

complies  with  the  upper  boundary  condition  (iii)  (a)  of 

§  7,  i.e.  (x)<  0,  K^,(x)>  0.  as  =  0,  it  follows 

o  0 

from  (6.16)  that  ^21^x^  =  0  and  from  (6.12)  that 

(  ^ 

>Vo,l30  -  Au(J)e^.ft  Kt(>)  c(u  03-'2) 

L  ? 

Hence 

X„(Oj  =  i  K,lO)  CM3) 


and  (13.9)  becomes 


1 


-  t>v  - 


[?(*,*  )J''  =■  I  -  [i-i.K,l0)]*/Hl<3) 

-  1  -  U  +  i«.*) 


by  (12.23). 


For  an  atmosphere  of  constant  scale  height  H,  = 

Gxx  =  0  by  (12.24),  and  (13*11)  and  (13.14)  give 


(13  *t) 


Fig.  lu  shows  (13.1b)  and  (13.16)  plotted  on  an  Argand 
diagram  lor  a  range  oi'  values  of  h/h.  The  only  physical 
quantity  that  enters  into  the  calculation  of  this  curve  is 
y/,  the  ratio  of  specific  neats  ior  air,  which  is  taKen  to 
be  7/b.  The  semicircular  part  of  the  plot  has  radius  K. 
and  L  is  the  point  corresponding  to  the  leading  lunar 
semi-diurnal  mode  with  H  =  7.6  km.  At  the  origin 


ti  -  1  =  (d  -  «il!  Y  (13.1 7) 

L  a  ^  k  I 


by  (13.16).  By  (3*b),  equation  (13.17)  gives 


ft  =  V  H 


(13-lS) 


I 


Fig. 


10  Argand  diagram  of  q  for  a  range  of  values 
of  h/H,  where  scale  height  H  is  constant. 

Y  =>  7/5  and  4  <■  =4(  ^  -1 )/  V  *  8/7. 


M 


against  h/H,  where  1J  is  shown  in  Fig.  10 


Fig.  11 
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Fig.  11  shows  the  magnification  (r>|  as  a  function  of  h/h. 

Trie  main  feature  of  tnis  curve  is  the  free  oscillation  at 
h/H  =  V  :  for  H  =  7*6  km  the  free  oscillation  arises  when 
h  =  10. b4  km. 

For  a  realistic  scale  height  profile  (fig.  I )  the 

evaluation  of  r>  by  (15. o)  or  (15.  a ;  generally  requires 

integrations  of  trie  homogeneous  equation  (5.1)  over  a 

range  (0,  x  .)  for  a  value  x,,  such  that  WKne  solutions  are 
111 

valid  for  x  >•  x^..  For  the  scale  height  profile  of  -Hig.  1, 
t]-1  is  snown  in  Fig.  12  lor  a  range  of  values  of  h  in  urn . 
(The  curve  between  h  *  7*5  and  9  km  is  omitted  to  avoid 
confusion  with  the  part  of  tr;e  curve  between  h  =  4  ana 
o  km  which  it  follows  closely).  On  comparing  Figs.  10 
ana  12  (or  Figs,  li  and  15),  tne  change  from  a  constant  to 
a  realistic  scale  height  is  seen  to  drastically  effect  the 
response  for  values  of  h  from  about  4  to  d  km  including 
that  of  the  leading  lunar  semi-diurnal  mode  at  L  (Figs. 

10  ana  12). 

In  Fig.  15  the  magnification  \{]\  has  a  very  sharp  peax 
at  h  =  6.77  km  with  a  magnification  of  16.5.  Tne  larger 
peak  at  h  =  10.5  km  is  similar  to  that  for  a  constant  scale 
heignt  (Fig.  11)  but  the  peak  magnification,  although  very 
large,  is  no  longer  infinite.  The  two  peaxs  of  Fig.  15 
nave  featured  in  previously  puDiished  response  curves 
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(wee  Kes  &  rtilxes ,  194/;  Jaccnia  &  Kopal ,  l'.-92 ;  J  iwa  ,  x  -o’.'; 
Hollingsworth,  1971).  uiwa  (iyt>d)  obtair.ea  a  maxima.:  at 
h  =  .•5.4  km  witn.  a  magnif icat ion  of  1.7;  and  a  third  maximum 
is  obtained  in  the  present  calculations  at  h  =  5.1  kx,  tut 
it  is  too  small  to  be  apparent  in  rig.  1'5  although  1 t  is 
shown  b.y  the  tight  loop  in  Fig.  12.  A  still  smaller 
maximum  obtained  by  Giwa  (lyoa)  at  h  =  1.9  ks  wouia  appear 
to  correspond  to  tne  ini’ 1  exion  that  is  apparent  in  the  curve- 
in  Fig.  12. 

As  multiple  reflexions  do  not  arise  with  a  constant 
scale  height  a  more  realistic  result  for  might  be 

expected  to  be  shown  by  Fig.  10  tnan  Fig.  12,  in  spite  of 
tr.e  scale  height  approximation  involved.  Using  tne  equation 
of  L.ie  locus  snown  in  Fig.  10,  Chapman  &  Lindzen  (1970) 
pointed  out  that  both  the  observed  lunar  semi-diurnal 
pressure  amplitude  of  ^  V’O  p-  b  and  phase  of  «•-  72°  are 
quite  well  represented  b.y  taking  H  =  6.7b  km.  The  predicted 
phase  can  be  read  off  directly  from  Fig.  10  as  the  an-*le  between 
OL  and  the  downward  axis,  where  0  is  the  origin.  L  is  plotted 
in  Fig.  10  for  H  *  7.6  km  and  a  phase  of  6d°  is  then  obtained. 

Associated  with  an  externally  applied  tiua :  potential 
Jlot.  is  an  additional  potential  k  A  arising  from  deformation 
of  the  Farth  and  a  vertical  velocity  component  at  tr.e  earth's 
solid  surface  of  -n'J?Q£/go,  where  h'  and  k  are  bove ' numbers. 
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For  ocean  regions  an  additional  vertical  velocity  r 


* i  ,,  is  introduced  by  the  ocean  tide  and  there  is  also  a:, 
additional  potential  .fl  ^  due  to  deiormation  ol  me  ocean 
(Hollingsworth,  1971).  11'  the  individual  excitations  are 

taken  to  be  additive  we  have 

Jl0  ~  -t  ftoc 


W0  -  WoC  -  £'  JlcS  1 2c 


(>3. 


Then  for  a  single  mode  by  (2.21)  and  (2.22) 

SI  =  0*0  Jlr  i-  jlc 
yi  =  Wc  t  i  %  J1E 


Cis.aO 


II  the  surface  oscillations  generated  by  (15.20)  and  ( 1^.21 ) 
are  denoted  by  PqJX  ,  PqW »  the  oscillation  recorded  by  a 
iand-based  barometer  is 

Z  =  *  ToW  +  (v/V'"’ 

Tne  third  term  on  the  right-hand  side  is  tne  pressure 
oscillation  due  to  the  small  vertical  displacement  of 
the  barometer  -h'/L,yg  .  For  a  single  mode  il5. tecores 
o y  (2.21)  and  (2.22) 


P  =  \  ♦  Tw  T  \  J1B  ae/Hl®\ 


('5.  2%) 


Hence  by  (1 0.2),  (10.12)  and  (10.15)  witn  x.  =  0 

/i 


?-  L.d0_?  Ve 

i  Hlo) 


^|NW  iO(o)[(u-fe--R,)jle 

f  N^IO  v)J}(u)duj 


V  IU4-) 


If  the  contributions  of  the  oceans  and  atmospheric  dissipation 
are  neglected  we  nave 


T  -  -  (K£-r)j\(\lwU>oK/£  (.  5  15) 


(15.25)  shows  that  the  elasticity  of  the  Earth  nas  the  effect 
of  multiplying  and  hence  the  resulting  tidal  pressure 
amplitude  by  (1+k-h'),  i.e.  0 y  about  0.70,  whicn  is  the  result 
i';i  veil  by  Hollingsworth  (1571). 


14.  Discussion 

Atmospheric  tidal  theory  has  been  aeveiopeu  and  applied  ov 
many  decades,  but  a  comprehensive  analytical  treatment  aion.: 
the  lines  of  that  of  butler  5c  Email  ( 1 5b  3 ) ,  which  essentially 
involves  the  forming  of  Ji’een's  functions,  has  been  lacking. 

A  difficulty  that  soon  arises  in  any  such  analysis  is  that  of 
manipulating  rather  cumbersome  equations.  Close  attention 


nas  therefore  neen  given  here  to  the  choice  of  notation  and 
formulation:  in  particular  it  seemed  wortrwr.ile  lntroducinr 


the  symbol  A  lor  the  cross-product  that  frequently  occurs. 

A  feature  of  the  analysis  has  seen  the  systematic  formulation 

of  results  for  different  atmospheric  variables,  i.e.  wind 

components,  temperature  and  pressure.  s regress  in  this 

direction  was  helped  by  tne  introduction  of  the  notation 

(Equ.  5.3)  and  the  derivation  of  15*10) •  another  feature 

of  the  analysis  has  been  tne  retention  of  the  arbitrary 
m _ 

constants  .T,  ,  ,  X'  until  they  could  oe  simultaneously 
o  1  oh 

eliminated  oy  the  introduction  of  two  boundary  conditions  to 
give  the  determinant  (8.1)  and  thence  the  general  relations 
of  §  8.  The  upper  boundary  condition  follows  earlier 
accounts  by  adopting  the  radiation  condition,  but  its 
application  in  terms  of  WKBJ  solutions  is  less  restrictive 
as  tne  actual  structure  of  the  atmospheric  'top'  does  not 
need  to  be  specified.  At  tne  .Lower  boundary  the  vertical 
component  of  velocity  is  retained  as  tne  usual  assumption 
of  a  zero  value  holds  only  for  a  surface-air  interface 
that  is  rigid  and  horizontal.  Undulations  o.'  t*  errain 
interact  with  the  primary  atmospheric  oscillations  and  set 
up  new  oscillations  for  whicn  vertical  velocity  components 
at  tne  surface  are  not  necessari i,y  zero.  ourface  tidal 
motions  are  another  source  of  excitation  at  the  lower 
uoundary  to  wnich  the  present  results  are  arplicable  in  terms 


of  vertical  motion 
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Application  01'  the  derived  formulae  to  quantitatively 
defined  sources  of  excitation  has  not  been  undertaken  within 
tne  scope  of  this  paper.  The  analysis  has  howe\Ter  teen 
developed  for  particular  heights:  the  first  case  at  the  lower 
boundary  (§  10)  leads  to  the  evaluation  of  surface  oscillation 
weighting  functions  for  a  range  of  solar  diurnal  and 
semi-diurnal  modes  (Equ.  10.  ib,  Figs,  d  to  •+).  The  second 
case  is  for  heights  above  the  region  of  excitation  ana  leads 
to  the  evaluation  of  tne  thermal  response  weighting  function 
V»t  for  the  same  modes  (Equ.  Ll.b,  Figs.  S  to  7),  rt  ,  W 
have  previously  been  presented  graphically  for  a  selection 
of  modes  and  provide  a  useful  means  of  understanding  the 
relationship  between  the  vertical  structure  of  a  thermal 
excitation  and  its  resulting  atmospheric  response  (troves, 
1975,  1976,  1977). 

Wo  attempt  nus  been  mace  to  define  the  \ imitations  of 
classical  theory  in  its  application  to  the  real  ntmospnere , 
but  reference  is  made  in  9  15  to  the  case  of  excitation  by 
the  leading  lunar  semi-diurnal  mode  of  tidal  potential  for 
which  classical  results  may  oe  unrealistic.  The  treatment 
of  the  surface  pressure  oscillation  due  to  a  tidal  potential 
(8  15)  has  been  of  long-standing  interest  and  provides  an 
example  of  the  application  of  (10.15):  tne  results  obtained 
are  gratifyingly  in  close  agreement  with  previous  accounts. 
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